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ABSTRACT 

Elliptical galaxy nuclei are the sites of the largest black holes known, but typ- 
ically show little or no nuclear activity. We investigate this extreme quiescence 
using Chandra X-ray Observatory observations of the giant elliptical galaxies 
NGC 1399, NGC 4472, and NGC 4636. The unique Chandra imaging power 
enables us to place upper limits of 7.3, 15 and 28 x 10^^ the Eddington lumi- 
nosity for the ~ 10^ - 10^ M© black holes in NGC 1399, NGC 4472, and NGC 
4636, respectively. The corresponding radiative efficiencies in this band are 4.1, 
24, and 620 x 10~^ using Bondi accretion rates derived from the Chandra hot 
interstellar gas surface brightness profiles. These limits are inconsistent with ba- 
sic advection-dominated accretion fiow (ADAF) models for NGC 1399 and NGC 
4472, indicating accretion onto the black hole at < 10% of the Bondi rate. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: 
elliptical and lenticular — galaxies: nuclei — X-rays: galaxies 
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1. Introduction 

Knowledge about the frequency and masses of supermassive black holes (SMBH) in 
galactic nuclei, and the characteristics of their associated radiation, deepen our understand- 
ing of galaxy formation and clarify its connection with the active galactic nuclei (AGN) 
phenomenon (e.g., Blandford 1999). While the intense energetics and rapid variabihty in the 
most luminous active galaxies form the empirical foundation of the SMBH/ AGN paradigm, 
such extreme objects are rare - especially in the local universe. 

Nevertheless, there is now a consensus that most, if not all, galaxies host SMBH. This 
is based, in part, on the high frequency of low-level nuclear activity in nearby galaxies (Ho, 
Filipenko, & Sargent 1997), and on statistical arguments connected to the cosmic X-ray 
background (e.g., Fabian & Iwasawa 1999). However, the strongest evidence for the ubiquity 
of SMBH comes from dynamical studies. Magorrian et al. (1998) found that, of 32 galaxies 
that could be fit with an axisymmetric dynamical model, 26 (29) require the presence of a 
central massive (~ 10^ — 10"'^'^ Mq) dark object at 95% (68%) confidence. More unambiguous 
and accurate measurements based on Hubble Space Telescope (HST) stellar or ionized gas 
spectroscopy confirms the presence of SMBH while deriving systematically lower SMBH 
masses (Gebhardt et al. 2000a; Merritt & Ferrarese 2001a). 

Although kinematically-studied luminous AGN are found to host SMBH (Gebhardt et 
al. 2000b; McLure & Dunlop 2000), the converse is not true: activity in most SMBH host 
galaxies is modest or absent. This can be quantified by considering the ratio of the nuclear 
luminosity to the Eddington luminosity (I/edd — 1-25 IO^^Msmbh, where the black hole mass 
-^SMBH is measured in solar units). The most powerful AGN have ratios on the order of 
unity, as expected from the theory of thin disks accreting at their maximal rate. However, 
this ratio is less than 10~^ and 10~^, respectively, in the two most securely measured SMBH 
- SgrA* at the Galactic Center (Quataert, Narayan, & Reid 1999) and NGC 4258 (Gammie, 
Narayan, & Blandford 1999). If these SMBH accrete at the Bondi (1952) rate, they must do 
so in a fundamentally different - and much less radiatively efficient - way than in luminous 
AGN. The theory of advection-dominated accretion fiows (ADAF; e.g., Narayan & Yi 1994), 
provides a predictive formalism for explaining the observed quiescence in the nuclei of SgrA*, 
NGC 4258, and other systems (Fabian & Rees 1995; Mahadevan 1997). 

Elliptical galaxies are the sites of the most luminous quasars, and also host the most 
extreme cases of quiescent SMBH. Although they contain nuclear black holes with the highest 
estimated masses and are observed to contain an ample supply of fuel in the form of X-ray 
emitting hot interstellar gas, they often emit less than 10~^Ledd (Loewenstein et al. 1998; 
Di Matteo et al. 2000). Radio/X-ray observations of some systems are inconsistent with 
predictions of the simplest models of an ADAF accreting at the Bondi rate, but can be 
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accommodated by variations that include outflow (Blandford & Begelman 1999; Quataert & 
Narayan 1999; Di Matteo et al. 1999) or convection (Narayan, Igumenshchev, & Abramowicz 
2000; Quataert & Gruzinov 2000). 

Prior to the launch of the Chandra X-ray Observatory, X-ray searches for nuclear point 
sources in elliptical galaxies were complicated in the imaging domain by centrally concen- 
trated hot gas that dominates the soft X-ray flux; and, in the spectral domain by an ex- 
tended hard component consisting of X-ray binaries. Thus the previous sensitivity to SMBH- 
associated X-ray emission, for even the nearest giant ellipticals, was ~ lO''^ erg s^-"^ (Roberts 
& Warwick 2000; Sulkanen & Bregman 2001), as neither the ASCA SIS/GIS nor the ROSAT 
PSPC/HRI had the necessary combination of sensitivity and spatial and spectral resolution 
to probe more deeply. Allen, Di Matteo, & Fabian (2000) claim to detect a very hard nuclear 
X-ray component in six elliptical galaxies with characteristics consistent with ADAF/outflow 
models, given the observed radio emission (Di Matteo et al. 1999). However, these results 
are based on a non-unique decomposition of ASCA spectra. 

Because of the presence of bright stellar cores and dust, the best approach for con- 
straining the level of SMBH activity is via high angular resolution X-ray observations (e.g., 
George et al. 2001): for a typical AGN X-ray/optical ratio, the Chandra detection thresh- 
old corresponds to a B magnitude ~ 24, while elliptical galaxy central surface brightnesses 
range from ~ 11 — 16 mag arcsec"^. This is further demonstrated by deep Chandra surveys 
in combination with optical follow-up. Low- luminosity AGN, most not optically detected, 
dominate the X-ray background; and, nonthermal emission is found to be common - at any 
given time ~ 10% of all bright galaxies host Lx ^ 10^^ erg s^^ AGN (Mushotzky et al. 2000; 
Barger et al. 2001). There is also a high frequency of nuclear X-ray point sources associated 
with local AGN, though often at levels that overlap in luminosity with X-ray binaries (Ho 
et al. 2001). Chandra has detected unresolved nuclear emission associated with quiescent 
SMBH in M31 at ~ 5 x lO-^Ledd (Garcia et al. 2000), the Galactic Center at ~ 10-^^Ledd 
(Baganoff et al. 2001), and NGC 6166 at ~ 10-^Ledd (Di Matteo et al. 2001). 

In this paper we use Chandra imaging data to place luminosity limits on any nuclear 
point sources in the elliptical galaxies NGC 1399, NGC 4472, and NGC 4636. These limits 
raise new questions about the radiative efficiency and rate of accretion onto supermassive 
black holes. Following Merritt & Ferrarese (2001b), we adopt distances of 20.5, 16.7, and 
15.0 Mpc for NGC 1399, NGC 4472, and NGC 4636, respectively. 
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2. Data Analysis and Results 

2.1. Observations, Data Reduction and Analysis 

NGC 1399, NGC 4472, and NGC 4636 were observed with the Chandra ACIS-S detec- 
tor for 40-60 ksec during Cycle 1. Results on the diffuse hot gas emission and the X-ray 
binary population, and a more detailed data reduction description, are included elsewhere 
(e.g, Angelini, Loewenstein, & Mushotzky 2001). The latest revised standard pipeline pro- 
cessing was used for NGC 4472 and NGC 4636, while an earlier version in combination with 
additional bad pixel removal, gain map correction, and quantum efficiency map correction 
was used for NGC 1399. Final useful exposure times are 56, 40, and 52 ksec for NGC 1399, 
NGC 4472, and NGC 4636, respectively. NGC 1399 was approximately centered on the S3 
chip; NGC 4636 and NGC 4472 were positioned at the standard ACIS-S aimpoint. 1.25-7 
keV images of the central 50" x 50" are shown for all three galaxies in Figures la, lb, and 
Ic. The positions of the optical nuclei from archival HST NICMOS (NGC 4472, NGC 4636) 
or WFPC2 (NGC 1399) images are marked. The brightest point sources in the images are 
easily visible and indicated in the images; their counting rates are in the 0.005-0.008 cts s~^ 
range. The bright X-ray core in NGC 4472 is resolved and has a thermal spectrum. Only in 
the case of NGC 4636 is there a possible coincidence of the nucleus with a Chandra X-ray 
source (just to the southwest of the nuclear position indicated in Figure Ic), with a 0.3-7 
keV count rate ~ 0.004 cts s~^ (see below). 0.3-7.0 keV surface brightness profiles are shown 
in Figure 2. For illustration purposes, the NGC 4636 profile is centered on the near- nuclear 
source. No excess emission from unresolved nuclear emission is apparent for the other two 
galaxies. 

The position of the nucleus in NGC 1399 is identified based on registration of globular 
cluster X-ray sources with their WFPC2 counterparts (Angelini et al. 2001); the resulting 
HST and Chandra centroids are consistent. Similarly, the nucleus of NGC 4636 is located 
based on registration with a bright NICMOS source and with the USNO catalog (Monet 
et al. 1998), and is found not to be consistent with the position of the near- nuclear source 
(Figure Ic). For NGC 4472, we adopt the unregistered i/ST centroid that is consistent with 
the peak in the soft thermal emission seen in the broadband Chandra image (Figure lb). 
3a upper limits to any nuclear X-ray point source emission are derived following Gehrels 
(1986). Counts are accumulated in inner ("source") and outer ("background") concentric 
annuli with radii that are varied to obtain the tightest upper limit. The Chandra calibration 
database enclosed energy function at the appropriate off-axis angle is used, and an exposure 
map correction is apphed for NGC 1399, which was centered off-axis. The assumed position 
of the SMBH is allowed to vary in a 3" x 3" box (corresponding to a very conservative 
adopted Chandra/ HST relative systematic pointing uncertainty) centered on the assumed 
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nuclear positions and the highest upper hmit in this box is adopted. These boxes encompass 
the peaks of the diffuse emission for NGC 4472 and NGC 1399. The diffuse hot gas emission 
is imphcitly assumed to be locally fiat on average when treated as "background" , consistent 
with the flat central surface brightness profiles (Figure 2). Luminosity upper limits are 
derived both in the 0.3-7 keV and 2-10 keV energy bands assuming a slope 1.5 power-law 
spectrum - the limits increase by ~ 2 for slope 0.5. Due to reduced dilution from the 
diffuse hot gas emission, the 2-10 keV count rate limits are 3-5 times lower than in the 
softer band; however, because of the declining effective area the energy flux upper limits are 
comparable. The 0.3-7 keV detection limit varies from a few 10^^ erg s~^ in regions of low 
surface brightness hot ISM emission to a few 10^^ erg in the core of NGC 1399. 

The 2-10 keV limits - all derived within 1" of the optical nucleus - are displayed in 
Table 1; these correspond to 28, 23, and 16 counts for NGC 1399, NGC 4472, and NGC 
4636, respectively - comparable to the brightest sources in the subimages shown in Figure 
1. Table 1 also includes less conservative, but perhaps more realistic, limits calculated by 
fixing the positions at the assumed nuclear positions described above. 

The inner gas density profile is estimated by fitting the point-source-subtracted surface 
brightness profiles with /3-models, 



out to some maximum radius, followed by analytical deprojcction with conversions obtained 
from spectral fits to the inner regions (best-fit temperatures kT = 0.8, 0.8, and 0.6 keV 
for NGC 1399, NGC 4472, and NGC 4636, respectively). The quality of fit dechnes as 
the maximum radius is increased due to deviations from the /3-model, but the profiles are 
sufficiently concentrated that the central density estimates - 4.6, 3.2, and 1.1 x 10~^^ gm 
cm-=^ for NGC 1399, NGC 4472, and NGC 4636, respectively - are robust. 



Individual black hole masses for the three galaxies studied here have been estimated by 
Magorrian et al. (1998) from two-integral dynamical modeling of ground-based kinematical 
data. A more accurate value can be derived using the correlation with bulge velocity dis- 
persion, (7, derived for galaxies with limits based on iJS'T spectroscopy. Therefore, we adopt 
SMBH masses from Merritt & Ferrarese (2001b) (Table 1); these are 3-5 times lower than 
in Magorrian et al. (1998). The resulting 2-10 keV luminosities in units of the Eddington 
luminosity and the Bondi (1952) accretion rates for adiabatic spherical infiow are displayed 




(1) 



2.2. Accretion Rates and Eddington Ratio Limits 
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in Table 1. Our derived upper limits on Lx/Ledd are 7.3 x 10~^, 1.5 x 10~^, and 2.8 x 10~^ 
for NGC 1399, NGC 4472, and NGC 4636, respectively. The corresponding efficiencies, 
assuming accretion at the Bondi rate Msondi, are 4.1 x 10~^, 2.4 x 10~^, and 6.2 x 10~^. 



3. The Nature of the Accretion Flow 

The conservative upper limits we derive are 1-2 orders of magnitude lower than the 
luminosities estimated by Allen et al. (2000) using ASCA spectral decomposition. Direct 
Chandra images (Angclini ct al. 2001) indicate that most of the hard X-ray emission spatially 
unresolved by ASCA originates in X-ray binaries. Since upper limits on the nuclear X-ray 
emission and estimates of the Bondi accretion rate have each decreased by ~ 10 — 30 from Di 
Matteo et al. (2000) (the latter because of lower black hole mass estimates)^, a re-examination 
of the nature of the accretion flow is in order. 

For our estimated black hole masses and Bondi accretion rates, we have computed 
spectra using the techniques described in Quataert & Narayan (1999); we aim to test the 
basic ADAF model and so do not include outflows or convection. The resulting predicted 2- 
10 keV luminosities, Ladaf, are^2x 10^\ 10^°, and lO^^ ergs s"^ for NGC 1399, NGC 4472, 
and NGC 4636, respectively (Table 1). Our observational upper limits are thus significantly 
below the ADAF predictions for NGC 1399 and NGC 4472, but consistent for NGC 4636. 
The predicted X-ray luminosities in the standard ADAF models for the cases of NGC 1399 
and NGC 4472 are dominated by Compton scattered synchrotron emission that originates 
close to the black hole. The luminosity is primarily sensitive to the gas density (i.e., accretion 
rate) and electron temperature. A factor of ~ 10 decrease in the accretion rate below the 
Bondi value can reduce the predicted luminosity to the observational limits. 

With the revised black hole masses and accretion rates we flnd that the radio flux 
predicted by the ADAF model is less than the observed flux for NGC 1399 and NGC 4636, 
but is in excess for NGC 4472; thus the radio flux problem highlighted by Di Matteo et al. 
(1999) is somewhat ameliorated. Finally, a preliminary HST measurement of the nuclear 
FUV emission in NGC 1399 (O'Connell et al., in preparation), - i/Li, ~ 7 10'^^ erg s^^ - also 
conflicts with pure-inflow ADAF models, although dust may be an issue at this wavelength. 

Our observations require a modiflcation of the standard ADAF/Bondi model such that 



•^For NGC 4636 our Bondi accretion rate is ~ 3000 times smaller than that of Di Matteo et al. (2000). 
A factor of ~ 10 is due to our smaller black hole mass; the remainder is primarily due to a lower inferred 
central gas density that may be underestimated since the accretion radius is unresolved in this case. 
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the gas density near the SMBH is reduced by at least a factor of ~ 10. This could be due 
to the role of strong outflows or convection in an ADAF. Alternatively, the central regions 
of these galaxies may not be in steady state due to heating by supernovae or by episodic 
accretion onto the central SMBH (e.g., Ciotti & Ostriker 1997). In addition the cooling time 
at the accretion radius -RBondi (Table 1) is < 10^ yr for all three galaxies; so, inflow with 
mass deposition could be important (Brighenti & Mathews 1999). 

4. Conclusions 

We have used Chandra observations to place new and unprecedented upper limits on the 
X-ray emission associated with ~ 10*^ — 10^ Mq nuclear black holes in three giant elliptical 
galaxies. The black holes clearly arc not fuel-starved: all have central hot gas densities 
exceeding 10~^^ gm cm^'^. Adopting masses from the Msmbh — cr correlation, and assuming 
that mass is inflowing at the Bondi (1952) rate, we derive 3a 2-10 keV luminosity upper limits 
of 7.3, 15 and 28 x 10~^LEdd and radiative efficiencies in this band of 4.1, 24, and 620 x 10~^ 
for NGC 1399, NGC 4472, and NGC 4636, respectively. Unless the black hole masses 
are considerably lower than predicted by the Msmbh — cr correlation, these observations 
are inconsistent with the simplest ADAF models in NGC 1399 and NGC 4472, and imply 
accretion onto the SMBH at below the Bondi rate. 

If the decline in the accretion rate is due to outflow, this raises the question of the 
ultimate fate of this gas as it encounters the dense hot interstellar medium in the core. Many 
questions about stability, intcrmittency, and regulation of gas flows in the inner regions of 
ellipticals remain to be addressed. The Chandra isophotes in the cores of NGC 1399 and 
NGC 4472 appear round and undisturbed, although there is an apparent spiral-like structure 
in the NGC 4636 diffuse emission. 

The more stringent limits required to further constrain the nature of the central accretion 
flow will be difficult. Our observations indicate > 0.01 super-Eddington X-ray binaries 
per squarc-arcsccond near the nucleus; so that detections below ~ 10^^ erg s~^ cannot be 
unambiguously identified as AGN. Although the hot gas emissivity precipitously declines at 
energies above 2 keV, so does the Chandra sensitivity, restricting the detection limit in gas- 
rich systems. Tighter hmits can be obtained on gas poor systems that, although of interest, 
may be fuel-starved and would thus shed little light on the accretion process. It would be 
especially useful to spectrally decompose the inner ~ 1" to place a limit on the fraction 
of non-ISM emission. This would appear as a hard excess and would require a very deep 
exposure due to the declining Chandra effective area mentioned above. 
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Table 1. Galaxy and Accretion Flow Characteristics 



galaxy 


d 


-^SMBH 


^Bondi 


-^Bondi 


-^edd 


T a 
-^Bondi 


7" b 
-t^ADAF 




T d 


NGC 


Mpc 


10« Mq 


II 


Mo yr-1 


erg s"^ 


erg 


erg s"^ 


erg s"^ 


erg 


1399 


20.5 


10.6 


0.36 


4.0 10-2 


1.3 10^^ 


2.3 


2 10^^ 


< 9.7 lO^s 


< 9.7 10^8 


4472 


16.7 


5.65 


0.24 


7.9 10-3 


7.1 lO^'' 


4.5 10^3 


1040 


< 6.4 10^8 


< 4.9 10^8 


4636 


15.0 


0.791 


0.049 


8.0 10-5 


9.9 10^5 


4.5 10^^ 


1036 


< 2.7 10^8 


< 1.8 10^8 



-O.lMBondiC^ 

''approximate expectation of standard ADAF model; see text 
^^2-10 keV upper limit from this paper for 3" box 
^^2-10 keV upper limit at assumed optical nucleus 
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Fig. 1.— Chandra ACIS-S3 1.25-7 keV images of the central 50" x 50" of (a) NGC 1399, (b) 
NGC 4472, and (c) NGC 4636. The open, 1" x 1" box indicates the position of the optical 
(HST) nucleus, the letter "A" the brightest point source. Relative shifts based on image 
registration have been applied for NGC 1399 (1.5")and NGC 4636 (0.9"). 
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Fig. 2. — Chandra ACIS-S3 0.3-7 keV central diffuse emission surface brightness distributions 
in cts arcmin"^ s~^. Annuli are centered on the X-ray centroids for NGC 1399 and NGC 
4472, and on the near nuclear source for NGC 4636. 



